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Abstract 

We present a supersymmetric grand unification model based on SO(IO) group with 5*4 flavor 



symmetry. In this model, the fermion masses are from Yukawa couplings involving 10 and 126 
Higgs multiplets and the flavor structures of mass matrices of both quarks and leptons are deter- 
mined by spontaneously broken 5*4. This model fits all of the masses and mixing angles of the 
quarks and leptons. For the most general CP-violation scenario, this model gives sin^is a wide 
range of values from zero to the current bound with the most probable values 0.02 — 0.09. With 
certain assumptions where leptonic phases have same CP-violation source as CKM phase, one gets 
a narrower range 0.03 — 0.09 for sin^ia with the most probable values 0.04 — 0.08. This model gives 
leptonic Dirac CP phase the most probable values 2 — 4 radians in the general CP-violation case. 
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I. INTRODUCTION 



The discovery of no- zero neutrino masses and lepton mixings have raised hoye to under- 
stand the mystery of flavor structures of quarks and leptons in a unified way[l]. Although 
many similarities between leptons and quarks make such unification plausible, the mixing 
pattern in the lepton sector is very different from that in the quark sector. However, there 
are now many grand unification models based on 50(10) gauge group that can give small 
quark mixings and large lepton mixings along with all their masses with few assumptions 

Another interesting possibility is that there may exist horizontal underlying flavor sym- 
metry. This is favored by leptonic mixing pattern with near maximal atmospherical mixing 
angle and the vanishing ^^13. A permutation symmetry between /i neutrino and r neutrino 
in the flavor basis has been proposed in recent years . Even though there is no apparent 
evidence of such symmetry in charged lepton and quark sector, it has been shown that the 
unified description of quarks and leptons with this symmetry is possible|4]. Applications of 
higher permutation group 5*3, 5*4 and AA to flavor symmetry also have been discussed in 
the literature 5l [Gil? I . Other discrete groups such as dihedral group DA and D5 have been 
studied |8i|9|. 

In this paper, we focus on the group SA x SO{10). SA has certain good features to be a 
flavor symmetry. First, it has three dimensional irreducible representation to accommodate 
the three generations of fermions naturally. Note that this is different from S3 because 
the largest irreducible representation of S3 has dimension two and therefore we have to 
treat one family of fermions different from other two. Second, it can be embedded into 
continuous group SU{3) or 50(3) [10]. As we will show below, SA symmetry also gives 
degenerate spectrum of the right-handed neutrinos naturally, which has some interesting 
consequences for the neutrino phenomenology. For example, in this case, one can use the 
resonant enhancement of leptogenesis for (quasi-)degenerate right-handed neutrinos to gen- 
erate enough baryon asymmetry ll|] . With the degenerate heavy right-handed neutrinos, 
the low energy neutrino flavor structure is determined by Dirac mass matrix at the seesaw 
scale completely, which makes it easier to reconstruct high energy physics from low energy 



observables. Some work has been done in this direction. In Ref . [12] , Lee and Mohapatra con- 
structed a 5*4 X 50(10) model, which naturally gives quasi-degenerate spectrum of neutrinos 
masses with small solar angle, which already has been ruled out by large mixing angle MSW 
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solution to the solar neutrino problem. In principle radioactive corrections may amplify 
the solar angle and keep the other two angles unchanged, but generally this needs extreme 
fine-tuning of parameters at the seesaw scale to realize it. On the other hand, in a recent 



paper 



10| by Hagedorn, hinder, and Mohapatra, a low energy scale non-supersymmetric 



model is presented based on 5*4 fiavor symmetry, which can accommodate current neutrino 

into a SUSY GUT framework 
12|. In this letter, we address 



data. Our goal is to see if we can embed the model of Ref . [10 
without running into the small solar angle problem of Ref. 
this question and find that we can build a realistic model based on 5*4 x 5*0(10) with the 
proper choice of the parameter space. 

In this model, all the quarks and leptons of one generation are unified into a 16 spinor 
representation of 5*0(10) and the Yukawa coupling structures of three generations are de- 



termined by 54. We use 10 and 126 representations of 5*0(10) for Yukawa couphngs to 
account for all the fermions masses and mixing angles [3] l^- Even though in the most gen- 
eral CP-violation case this model has 18 complex parameters, it is not obvious whether it can 
accommodate all observed masses and mixing angles because of constraints from 54 fiavor 
symmetry and the correlations between quarks and leptons indicated by 5*0(10) unifica- 
tion. For instance, with the particle assignment of 54 in this model, the heavy right-handed 
neutrino mass matrix is proportional to an identity matrix, and the Dirac mass matrix of 
neutrino determines the mixing among light neutrinos completely. The general mechanism 
to generate the lepton sector mixing independently from the quark sector by right-handed 
neutrinos does not work in this model. On the other hand, one may argue that since the 
total number of parameters is much larger than that of obervables, this model may lose 
predicability even if it can fit all the obervables. We find this not to be the case. It turns 
out that half of complex phases can be rotated away by choices of basis and redefinitions 
of the right-handed fields of charged leptons and down-type quarks. For the most general 
CP-violation case, this model gives wide range of sin6'i3 from zero to current bound with 
the most probable values 0.02 — 0.09. The most probable values of leptonic CP phase are 
2 — 4 radians. With certain assumptions where the leptonic phases have same CP-violation 
source as CKM phase, one gets narrower predicted range 0.03 — 0.09 for sin ^13 with the 
most probable values 0.04 — 0.08. 

Some issues about Higgs sector still need to be addressed. As we have six 10s and three 
126s, without analyzing the 5*4 x 5*0(10) invariant Higgs potential, whether or not we can 
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get the desired vacuum configuration still remains an open question. We do not concern with 
doublet-doublet splitting and doublet-triplet splitting problems in this paper. With such 
rich Higgs fields, we assume they can be realized in some way. And another fact we should be 
careful is that generally the discrete flavor symmetry can enhance the accidental global sym- 
metry of Higgs potential and lead to unwanted massless Nambu-Goldstone bosons. There 
are ways found in the literature to avoid it. One can introduce gauge singlet Higgs flelds 
whose couplings are invariant under discrete symmetry but break the global symmetrv[l5|. 
or introduce soft terms which break discrete symmetry and global symmetry 161]. 

The paper is organized as follows: in Section 2, we present an SO (10) model with 5*4 
flavor symmetry and present the mass matrices of quarks and leptons; in Section 3, we 
present a detailed numerical analysis including CP violation in quark and lepton sector. We 
end with conclusions and remarks in Section 4. 



II. SUSY SO(IO) MODEL WITH S4 FLAVOR SYMMETRY 

The group 5*4 is the permutation group of the four distinct objects, which has 24 distinct 
elements. It has flve conjugate classes and contains flve irreducible representations 1,1', 2, 3 
and 3'. Our assignment of fermions and Higgs multiplets to 5*4 x 5*0(10) are shown in Table 



Fermions 


Higges Bosons 






Ao 


Al,2 


Ho 




-^3,4,5 


{3'} X {16} 


{1} X {210} 


{1} X {126} 


{2} X {126} 


{1} X {10} 


{2} X {10} 


{3} X {10} 



TABLE I: Transformation property of fermions and Higgs multiplets under 54 x 5'O(10) 

In this model, we assign three generations of 16 to 3' irreducible representation of S4, 
because 3' can be identifled with the fundamental representation of continuous group SO (3) 
or SU(3)[l^ {25!. In Higgs sector, because of 3' x 3' = 1 + 2 + 3 + 3', to make Yukawa coupling 
5*4 invariant, Higgs flelds can not belong to 1'. 1 is necessary for phenomenological reason, 
otherwise all of the mass matrices would be traceless. To get symmetric mass matrices which 
is required by group structure of 16 • 16 ■ 10 or 16 • 16 • 126, Higgs should not belong to 3'. 
We include both 2 and 3 to get realistic mass and mixing of quark and lepton. One might 
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think six 10 Higgs fields transforming as 1 + 2 + 3 under 5*4 are enougli. But tliere are two 
reasons wliy we also need 126, one is to give right-handed neutrinos heavy masses and the 
other is to fix the bad mass relation between quark sector and lepton sector indicated by 
16 ■ 16 ■ 10. In this sense, our choice of Higgs fields is minimal. 

The breaking of 5*0(10) to Standard Model(SM) can be realized in many ways. In this 
model, we choose 210 Higgs field, which is 1 under 5*4 transformation, to break SO{10) 
to SU{2)l X SU{2)n x SU{4:)c (0224) while keep the ^4 symmetry. We choose (1,3,10) 
components of only Aq (the numbers denote representation under the G224) to get vev vr 
that breaks G224 down to the SM and gives heavy masses to right-handed neutrinos. With 
this breaking pattern, 5*4 symmetry is kept down to the electroweak scale. 

To see what this model implies for fermion masses, let us first explain how the MSSM 
doublets emerge. Besides the SU{2)l Higgs doublets from submultimplets (2, 2, 1) and 
(2, 2, 15) contained in 10 and 126 respectively, we also have Higgs doublets contained 
in (2, 2, 10) © (2, 2, 10) from 210. Furthermore, to obtain anomaly-free theory, we need 
to introduce three 126, which we denote by A, that also contain Higgs doublets. Al- 
together, we have fourteen pairs of Higgs doublets: 0„ = (i7j„, Aj„, Aj„, $u2), 4>d = 
{Hid, ^jd, ^jd, ^di, ^d2), where i = 0, 5 and j = 0, ...2. As noted, six pairs from Hs, three 
pairs from As, three pairs from As and two pairs from $. We can write Higgs doublet mass 
matrix as (pu^H'P'd- can be diagonalized by XMhY'^ ■, which X and Y are unitarity 
matrices acting on (f)u and (pd respectively. At the GUT scale, by some doublet-triplet and 
doublet-doublet splitting mechanisms, we assume only one pair of linear combinations of 
X*^0u/3 and Y*p(l)di3, say X*^0„/3 and Yip(t)dfi, has masses of order of the weak scale and all 
others are kept super heavy near GUT scale, which generally can be realized by one fine- 
tuning of the parameters in the Higgs mass matrix. The MSSM Higgs doublets are given by 
this lightest pair: if^ssM ^ and if^^ssM ^ Yj*^0^^. Since we focus on the structures 

of Yukawa couplings, we do not discuss the details of the splitting mechanisms that lead to 
the above results. 

With Higgs fields and fermions listed in Tabl^ we can write down 5*4 x 50(10) invariant 
Yukawa coupling as [26] 

W^Yukawa = (^1^1 + ^2^2 + ^S^s) (/^O^O + /oAo) 
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+ _L(^2*2 - *3*3)(^1^1 + /2A1) + -^(-2*1*1 + *2*2 + *3^3)(/il^^2 + /2A2) 
+hs[{^2^3 + ^'3^2)^3 + (*1*3 + ^'3^1)^4 + (^'l^'2 + ^'2^'l)^5]- (1) 

After clcctrowcak symmetry breaking, (2,2,1) of Hi{i = 0, ...,5) component acquires 
vevs (denoted by (Hj)'^ and {Hi)'^). And (2, 2, 15) sub-multiplet of Aj(j = 0, 2) also get 
induced vevs. Their vevs are denoted by (Aj)" and (Aj)'^(j — 0, 1, 2). 

The mass matrices for the quarks and the leptons have following sum rules: 



Ml 



M. 



(10) 



3M 



(126) 



(2) 
(3) 
(4) 
(5) 
(6) 



where 



(10) 



(126) 



oo — 2a2 

as 
04 
60 - 262 

do - 2d2 



Co - 2e2 





05 

ao + oi + 02 
03 

bo + bi + 62 

&3 



(io + c^i + d2 




eo + ei + 62 




04 
03 

Oo — tti + a2 
h 

h 

bo-bi + 62 




do - di + d2 






(7) 



(8) 



(9) 



(10) 



eo - ei + 62 . 

and where and 6i are products of the type h{Hi)'^ and h{Hi)'^ respectively. Similarly, we 
use dj and Cj to denote products of the type /(Aj)" and fi^j)'^ respectively. The MSSM 
vevs are given by Vu = Xlp{(f)up) and va — Y^i^^cf)^^), where we use Vu and va to denote 
vevs of i/^ssM ^MSSM respectively. The Yukawa couplings and vevs of Higgs fields in 
general are complex, and there are 18 complex parameters. We choose a basis in which the 
down-quark mass matrix is diagonalized and set 63 = 0, 64 = 0, and 65 = 0. Note this is our 
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main difference witfi Ref. [12], wliere tliey clioose a basis in wliicli up-quark mass matrix is 
diagonal and set off-diagonal entries of to zeros, wfiicli leads to small solar mixing angle. 
In the basis we choose, the charged lepton mass matrix is also diagnolized. Therefore, the 
phases of 60, ^i, ^2, Cq, Ci, and ti can be rotated away by redefining 3 right-handed down- type 
quarks fields and three right-handed charged leptons. We treat 60, ^1, ^2, Cq, Ci, and 62 as real 
parameters in later analysis, and they can be determined by the masses of down-quark and 
charged lepton completely. 

Because the mass matrix of down-quark sector is diagnolized and M„ is symmetric, one 
can have 



CKM^^^uVCKM-, 



(11) 



where M„ = diag{mu,mc,'mf)). By fitting mass matrix of up-quark in Eq. fllll) . parameters 
as, a4, as can be determined. In addition, we get three conditions among the parameters 
ao,ai,a2,dQ,di, and d2- Therefore, there are three complex parameters left to be determined 
by masses and mixings of neutrino sector. Without loss of generality, we choose dQ,di, and 
(^2 to be determined by fitting of neutrino sector. And Dirac neutrino mass matrix can be 
written conveniently as 

'x 0' 



= VckmMuVckm - ^mt 



y 
z 



(12) 



with 



11 1 
X = — (rfo - 2(i2), y = — idQ + di + d2),z = — (do - rfi + ^2) 



(13) 



Because we know nothing about leptonic phases, in principle, there is no constraint on the 
phases of do,di, and d2. 



To see how this model can give a large atmospherical mixing a ng 



mate analysis first. Using first order Wolfenstein parameterization 
VcKM^uVcKM can be written as 

mt - A^X^{l-ir]- p) X^ + A^\^ ■■ 

AX'^il-ir]- p) -A\^ 1 



e, we give an approxi- 



17( 1 for the quark mixing. 



(14) 
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where we use rric/mt ~ A'* and rriu/mt ^ A*. Therefore, to get near maximal mixing of 623, 
y and z should satisfy 



A^(l + A) - % ^ 1 - 4^. 



(15) 



III. DETAILED NUMERICAL ANALYSIS 



To see if the model is phenomenologically acceptable, we first fit the masses of the charged 
leptons and down-type quarks using the mass values of leptons and quarks at the GUT scale 
with tan/5 = 10 [27] given in Ref.[18j: 



input observable 


tan/3 = 10 


ruu (MeV) 


0.7238t[i:lif 


rric (MeV) 


210.3273;^?;°^^ 


mt (GeV) 


00 /I Qoo+30.2676 
OZ.^OOO_i4 7g8g 


rud (MeV) 


i.OUJO_o.2304 


rris (MeV) 


29.9454tt^°°^ 


TTib (GeV) 




me (MeV) 


n Qe;oc:+0.0003 


(MeV) 


7^ 671 ^+'^■0^'^^ 


rrir (GeV) 


1 OQOO+0.0013 

i.zyzz_o.ooi2 



We use standard parametrization form for the Vckm and take the following values at 
I ' 

the scale Mz[l9i]: sin^gi2 = 0.2272, sin ^^13 = 0.00382, sin ^^23 = 0.04178 and the CP phase 
6q = ^, where we use subscript q to distinguish them from the lepton section mixing angles. 
And we use RGE running factor r] = 0.8853. At the GUT scale, we have the Vckm 



0.973841 0.227198 0.00169092 - 0.00292876i 

-0.227079 - 0.000134603i 0.97298 - 0.000031403i 0.0369876 
0.00675837 - 0.00284968z -0.0364044 - 0.000664834i 0.99912 



(16) 
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A. Quark and charged lepton sector 

Using the central values of charged lepton and down-quark masses at GUT scale, 
bo,bi,b2,eo, ei and 62 are solved from Eq.([S]) and Eq.Q (in Mev) 

60 = 387.756, 61 = -539.649, 62 = 193.27, Cq = -22.7734, d = 22.8717, 63 = -11.5298. 

(17) 

For up-quark sector, by solving Eq. lITT]) and Eq.(j2]), we get values of 03,04,05 and three 
conditions for oq, Oi, 02, do, di, d2 (in Mev): 

03 = -2990.72 - i54.757, 04 = 554.859 - i234.705, 05 = -66.748 + i8.155, 

oo - 2a2 + c/o - 2^2 = 14.628 - i3.162, oq + Oi + 02 + rfo + t^i + = 308.363 + i3.977, 

ao- ai + a2 + do- di + d2 = 82288.5 - i7.169 x 10"^ (18) 

We can see that accommodation of hierarchical structure of fermions masses is realized by 
adjusting the parameters, 5*4 flavor symmetry itself does not provide hints on it {2^. 



B. Neutrino sector 



In this model, the light neutrino mass matrix is given by type-I seesaw[21I]. The mass 
matrix of right-handed neutrinos is proportional to an identity matrix due to the 5*4 quan- 
tum number assignment, therefore the Dirac mass matrix determines the lepton sector 
mixing because the charged lepton mass matrix is diagnolized. 

1 



foVR 



J. V J. ^ J. V J. ^ . 



(19) 



This model gives hierarchical neutrino mass spectrum naturally. One can choose /o ~ 1 and 
vr ~ lO^^GeV, so the mass of the heaviest light neutrino is around 10~^ — lO^^eV. 

The fit of neutrino sector are found by scanning whole parameter space spanned by x, y 
and z under the constrain of the current experiment requirements. 

We choose the standard parametrization for the lepton sector mixing: 



U 



C12C13 



S12C13 



i5 



~C23Si2 — S23Si3Ci2e*'' C23C12 — S23Si3Si2e*'' 



■523C13 



■S23S12 ~ C23Si3Ci2e 



'■523^12 ~ C23Si3Si2e C23C13 



.diag(e-^^i/2^e-*^2/2^^) ^20) 
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with Cij = cos 6ij and Sij = sin 6ij. 5 is the Dirac phase and (fi, (p2 are Majorona phases 
of neutrinos. These phases have range from to 27r. 



We take 3a experiment bound j23j: 



0.24 < sin^^ia < 0.40 
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0.34 < sin^ 623 < 0.68 
sin^ ^13 < 0.040 
0.024 < Am|/Am^y^ < 0.040. (21) 

As mentioned earher x, y, and z generally are complex numbers. For the most general 
CP- violation case, we treat the phases of x, y, and z as random input numbers with range 
— 27r. The results are shown in Fig.(l). In this case, sin6'i3 has wide range from zero to 
the current bound with the most probable values 0.02 — 0.09 as shown in Fig. 1(a). Fig. 1(b) 
shows the correlation between sin ^23 and sin ^13. Fig. 1(c) is the value distribution of Dirac 
CP-violation phase in the lepton sector. The allowed range of 6 is quite large from to 27r 
radians with the most probable values 2 — 4 radians. Two Majorana phases ipi and (f2 have 
wide range from to 2tt as shown in Fig. 1(d), which is expected. 

Now we consider an interesting special case where x, y, and z are all real. Note the 
complexity of JoVr only contributes an overall phase to the light neutrino mass matrix, 
which can be rotated away. Therefore, in this case leptonic CP-violation phases have same 
source as CKM phase. 

The allowed range 0.03 — 0.09 for sin ^13 is narrower compared to the general case, and 
the most probable range is 0.04 — 0.08 as shown in Fig. 2(a). Unlike Fig. 1(b), Fig. 2(b) 
exhibits an interesting correlation between sin6'23 and sin6'i3. If we take the central value 
of ^^23 = f , we can get two much narrower ranges for sin6'i3. One is 0.055 — 0.06, and 
the other is 0.070 — 0.075. The values of 6 are 2.8 — 3 radians, and 6.0 — 6.1 with small 
possibility as shown in Fig.2(c). Fig. 2(d) shows the allowed values of two Majorana phases. 
Note this parameter region is just left-up corner of Fig. 1(d) for the most general case. The 
most probable value ranges for (fi and (f2 are 0.02 — 0.15 radians and 6.19 — 6.25 radians 
respectively. 

For illustration, we give a typical example of fit for this case. We take 

X = 0.0139726, y = 0.025914, z = 0.273173 (22) 
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25 





0.04 0. 

sin0i3 



(b) 




FIG. 1: Numerical analysis for the most general case where x, y, and z are complex consistent with 
current experimental bound Eq. ()2ip . (a) Value distribution of sin^ia. (b) Correlation between 
sin023 and sin^ia. (c) Value distribution of leptonic Dirac CP-violation phase, (d) Scatter plot of 
two Major ana CP-violation phases ^pi and ip2- 



and solve do, ^i, (^2, ^o, ai, a2 from Eq. f|T3|) and Eq. f|T8|) (in Mev) 

do = 8602.18, di = -10191.2, c/2 = 3725.19, Oq = 18935 + i0.271681, 

ai = -30798.9 + il.9887, 02 = 10036.1 + il. 71701. (23) 

With these parameters values as input, one then obtains for the neutrino parameters 

sin ^12 ~ 0.53, sin 6^23 ~ 0.73 

(24) 

sin ^13 ^ 0.054, Ami/ Am\^j^,j - 0-031- 

And light neutrino masses are rrii = 0.00774eV, 1712 = 0.01 18eV, = 0.051eV, which are 
normalized by Am'^i = 2.6 x 10^'^eV. The Dirac phase appearing in MNS matrix is 5 = 2.84 
radians. And two Majorona phases are (in radians): (fi = 0.093, !f2 = 6.21. The Jarlskog 
invariant [24| has the value Jcp = 1.80 x 10"'^. One can evaluate the effective neutrino mass 
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3 3.5 4 4.5 5 5.5 6 0.2 0.4 0.6 0.8 1 

5 Vl 



(c) (d) 
FIG. 2: Numerical analysis for case where x, y, z are real consistent with current experimental 
bound Eq. ipTI) . (a) Value distribution of sin 6*13. (b) Correlation between sin 6*23 and sin ^13. 
(c) Value distribution of leptonic Dirac CP-violation phase, (d) Scatter plot of two Majorana 
CP-violation phases ipi and 992 • 

for the neutrinoless double beta decays process to be 

\Y.Ulim^i\ ^ 0.009 eV. 

IV. SUMMARY AND CONCLUSION 

In summary, we build a supersymmetric 50(10) model with SA flavor symmetry. The 
three dimensional irreducible representation of 5*4 group unify three generations of fermions 
horizontally. 10 and 126 Higgs fields have been used to give the Yukawa couplings and 
generate all the masses and mixings of quarks and leptons. This model accommodates all 
obervables including CKM CP- Violation phase. We studied the prediction of this model 
in the neutrino sector. For the most general CP-violation case, this model gives the most 
probable values 0.02 — 0.09 for sin ^13. In a special case where leptonic phases have same 
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CP-violation source as CKM phase, one gets narrower range 0.03 — 0.09 for sin ^13 with the 
most probable values 0.04 — 0.08. 

le model we present here, the masses of light neutrinos purely come from the type-I 



In t 

seesaw 2l|]. Generally, one also can include the contribution from type-II seesaw[22], which 
can generate a scenario with degenerate neutrino mass spectrum naturally because of the 
5*4 symmetry if the type-II seesaw dominates the contribution to the light neutrinos masses. 
It is interesting to study the mixing pattern and its radioactive stability. We leave this 
possibility for future work. 

V. ACKNOWLEDGMENTS 

The authors thank R.N. Mohapatra for suggesting the problem and useful discussion, K. 
Hsieh, J.H. Kang and P. Rastogi for comments. H.B.Y. specially thanks C. Hagedorn for 
helpful comments on the discrete symmetry. This work is supported by the National Science 
Foundation Grant No. PHY-0354401. 



[1] For a recent review of the theoretical issues in neutrino physics, see R. N. Mohapatra and A. 



Y. Smirnov, arXiv:hep-ph/0603118, to appear in Ann. Rev. Nucl. Science, 56 (2006). 
[2] T. Fukuyama and H. Nishiura, |hep-ph /9702253[ R. N. Mohapatra and S. Nussinov, Phys. Rev. 
D 60, 013002 (1999); E. Ma and M. Raidal, Phys. Rev. Lett. 87, 011802 (2001); C. S. Lam, 



hep-ph/0104116i T. Kitabayashi and M. Yasue, Phys.Rev. D67 015006 (2003); W. Grimus 



and L. Lavoura, |hep-ph/0305046 ; 0309050; Y. Koide, Phys.Rev. D69, 093001 (2004);Y. H. 
Ahn, Sin Kyu Kang, C. S. Kim, Jake Lee, |hep-ph/0602160| A. Ghosal, |hep-ph/ 0304090; for 
examples of such theories, see W. Grimus and L. Lavoura, Phys. Lett. B 572, 189 (2003); 
W. Grimus and L. Lavoura, J. Phys. G 30, 73 (2004). 
[3] W. Grimus, A. S.Joshipura, S. Kaneko, L. Lavoura, H. Sawanaka, M. Tanimoto, 



hep-ph/0408123 ; R. N. Mohapatra, JHEP, 0410, 027 (2004); A. de Gouvea, Phys.Rev. D69, 
093007 (2004); R. N. Mohapatra and W. Rodejohann, Phys. Rev. D 72, 053001 (2005); T. 
Kitabayashi and M. Yasue, Phys. Lett,. B 621, 133 (2005); R. N. Mohapatra and S. Nasri, 
Phys. Rev. D 71, 033001 (2005); R. N. Mohapatra, S. Nasri and H. B. Yu, Phys. Lett. B 



13 



615, 231 (2005); R. N. Mohapatra, S. Nasri and H. B. Yu, Phys. Rev. D 72, 033007 (2005) 
Z. z. Xing, |arXiv:hep-ph/0605219l Z. z. Xing, H. Zhang and S. Zhou, |arXiv:hep-ph/0607091 



[4] K. Matsuda and H. Nishiura, Phys. Rev. D 73, 013008 (2006); A. Joshipura, [hep-ph/0512252 
R. N. Mohapatra, S. Nasri and H. B. Yu, Phys. Lett. B 636, 114 (2006). 

[5] L. J. Hall and H. Murayama, Phys. Rev. Lett. 75, 3985 (1995); C. D. Carone, L. J. Hall and 
H. Murayama, Phys. Rev. D 53, 6282 (1996); P. F. Harrison and W. G. Scott, Phys. Lett. 
B 557, 76 (2003); J. Kubo, A. Mondragon, M. Mondragon and E. Rodriguez- Jauregui, Prog. 
Theor. Phys. 109, 795 (2003) Erratum-ibid. 114, 287 (2005); S. L. Chen, M. Frigerio and 
E. Ma, Phys. Rev. D 70, 073008 (2004) Erratum-ibid. D 70, 079905 (2004); L. Lavoura and 
E. Ma, Mod. Phys. Lett. A 20, 1217 (2005); R. Dermisek and S. Raby, Phys. Lett. B 622, 
327 (2005); F. Caravaglios and S. Morisi, |arXiv:hep-ph/05¥3234 [; W. Grimus and L. Lavoura, 
JHEP 0508, 013 (2005); N. Haba and K. Yoshioka, Nucl. Phys. B 739, 254 (2006); S. Morisi, 



arXiv:hep-ph/0604106t R. N. Mohapatra, S. Nasri and H. B. Yu, Phys. Lett. B 639, 318 



(2006); R. Jora, S. Nasri and J. Schechter, hep-ph/0605069 S. Morisi, arXiv:hep-ph/0605167 



[6] S. Pakvasa and H. Sugawara, Phys. Lett. B 82, 105 (1979); E. Derman and H.-S. Tsao, Phys. 
Rev. D 20, 1207 (1979); D.-G. Lee and R. N. Mohapatra, Phys. Lett. B 329, 463 (1994); 
R. N. Mohapatra, M. K. Parida and G. Rajasekaran, Phys. Rev. D 69, 053007 (2004) ; E. Ma, 
Phys. Lett. B 632, 352 (2006); C. Hagedorn, M. Lindner and R. N. Mohapatra, JHEP 0606, 
042(2006). 

[7] E. Ma and G. Rajasekaran, Phys. Rev. D 64, 113012 (2001); K. S. Babu, E. Ma and 



J. W. F. Valle, Phys. Lett. B 552, 207 (2003); K. S. Babu and X. G. He, |arXiv:hep-ph/0507217 
A. Zee, Phys. Lett. B 630, 58 (2005); S. L. Chen, M. Frigerio and E. Ma, Nucl. Phys. B 724, 
423 (2005); G. Altarelh and F. Feruglio, Nucl. Phys. B 741, 215 (2006); X. G. He, Y. Y. Keum 
and R. R. Volkas, JHEP 0604, 039 (2006); B. Adhikary, B. Brahmachari, A. Ghosal, E. Ma 
and M. K. Parida, Phys. Lett. B 638, 345 (2006); E. Ma, H. Sawanaka and M. Tanimoto, 



arXiv:hep-ph/0606103t E. Ma, arXiv:hep-ph/0607190 



[8] W. Grimus and L. Lavoura, Phys. Lett. B 572, 189 (2003); W. Grimus, A. S. Jo- 
shipura, S. Kaneko, L. Lavoura and M. Tanimoto, JHEP 0407, 078 (2004); G. Seidl, 



arXiv:hep-ph/0301044 T. Kobayashi, S. Raby and R. J. Zhang, Nucl. Phys. B 704, 3 (2005). 



[9] E. Ma, arXiv:hep-ph/0409288 ; C. Hagedorn, M. Lindner and F. Plentinger,Phys. Rev. D 74, 



025007 (2006) 

14 



[10] C. Hagedorn, M. Lindner and R. N. Mohapatra, JHEP 0606, 042(2006). 

[11] J. R. Ellis, M. Raidal and T. Yanagida, Phys. Lett. B 546, 228 (2002); A. Pilaftsis and 

T. E. J. Underwood, Nucl. Phys. B 692, 303 (2004). 
[12] D. G. Lee and R. N. Mohapatra, Phys. Lett. B 329, 463 (1994) 
[13] K. S. Babu and R. N. Mohapatra, Phys. Rev. Lett. 70, 2845 (1993). 

[14] T. E. Clark, T. K. Kuo and N. Nakagawa, Phys. Lett. B 115, 26 (1982); C. S. Aulakh and 
R. N. Mohapatra, Phys. Rev. D 28, 217 (1983); C. S. Aulakh, B. Bajc, A. Melfo, G. Senjanovic 
and F. Vissani, Phys. Lett. B 588, 196 (2004); B. Bajc, A. Melfo, G. Senjanovic and F. Vissani, 
Phys. Rev. D 70, 035007 (2004); T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac and 
N. Okada, Phys. Rev. D 72, 051701 (2005); H. S. Goh, R. N. Mohapatra and S. Nasri, Phys. 
Rev. D 70, 075022 (2004). 

[15] T. Kobayashi, J. Kubo and H. Terao, Phys. Lett. B 568, 83 (2003). 

[16] K. Y. Choi, Y. Kajiyama, H. M. Lee and J. Kubo, Phys. Rev. D 70, 055004 (2004). 

[17] L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 

[18] C. R. Das and M. K. Parida, Eur. Phys. J. C 20, 121 (2001). 



[19] A. Hocker and Z. Ligeti, |arXiv:hep-ph/0605217 



[20] Q. Shafi and Z. Tavartkiladze, Phys. Lett. B 594, 177 (2004). 

[21] P. Minkowski, Phys. Lett. B67 , 421 (1977); M. Cell-Mann, P. Ramond, and R. Slansky, Su- 
pergravity (P. van Nieuwenhuizen et al. eds.), North Holland, Amsterdam, 1980, ; T. Yanagida, 
in Proceedings of the Workshop on the Unified Theory and the Baryon Number in the Universe 
(O. Sawada and A. Sugamoto, eds.), KEK, Tsukuba, Japan, 1979; S. L. Glashow, The future 
of elementary particle physics, in Proceedings of the 1979 Cargese Summer Institute on Quarks 
and Leptons (M. Levy et al. eds.). Plenum Press, New York, 1980, pp. 687; R. N. Mohapatra 
and C. Senjanovic, Phys. Rev. Lett. 44 912 (1980). 

[22] G. Lazarides, Q. Shafi and C. Wetterich, Nucl.Phys.B181, 287 (1981); R. N. Mohapatra and 
C. Senjanovic, Phys. Rev. D 23, 165 (1981). 

[23] M. Maltoni, T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 6, 122 (2004). See 



Appendix C of aiXiv hep-ph/0405172 for updated fit. 
[24] C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985). 

[25] If one gives up the possible embedding of 5*4 group to continuous group, one can choose 3 
and the mass matrices for fermions do not change. 



15 



[26] For the products and Clebsch-Gordan coefficients of 5*4 group, one can see for example the 
Appendix A of Ref.Q- 

[27] In this model, the value of tan/3 is not determined as in MSSM. We take tan/3 = 10 as an 
example. 

[28] One can use softly broken discrete flavor symmetry to understand fermion mass hierachy, see 



for example 



16 



